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Abstract 


Regarding to the conveyance of materials the assignment of stock-keeping units (SKUs) is one of the 
most important tasks in a warehouse. The units in the interest of their easier handling can be sorted 
into groups. To classifying SKUs the most frequently used method is the ABC analysis. The traditional 
approach implements a simplified scheme, which is mostly based on annual dollar usage. Classifying 
units according to one characteristic can lead to rushed result, because the assignment of stock-keeping 
units might be influenced by other factors, like number of orders, weight or lead time. In recent years 
a number of decision models have been developed in order to make it possible to consider more than 
one criterion in the same time. In general, in the case of representing a model the introduction of the 
algorithm is more preferable than the visualization of data. However the visualization of results can 
relieve the comprehension of the given problem and the way, how the model works, even in relation 
to sorting problems. The aim of this research is to introduce a plotting method which visualizes the 
results coming from the Case-based distance model developed by Chen, Kilgour and Hipel. This method 
strengthens the classification by providing separating nonlinear surfaces. 

Key words: case-based distance model, inventory management, visualization. 


Introduction 


Many obstacles can emerge during decision making and support. It can frequently occur 
that besides the large amount of data the decision has to be made by considering more than one 
factors or characteristics. For instance, the sorting problem in inventory management where 
thousands of stock-keeping-units have to be sorted into groups based on their characteristics. 
Concerning this problem numerous methods have been developed to simplify, analyze and 
relieve the multicriterial decision making situation 1n recent years. Visualization techniques can 
also help the decision support procedure by making the sort of problem more understandable. 

The traditional ABC analysis is the most frequently used approach to classify SKUs. The 
traditional approach implements a general, simplified scheme, which is mostly based on annual 
dollar usage. Besides its simplicity one of its advantages is the ease of representation of the 
results in the Pareto-chart (Harry, Mann, De Hodgins, Hulbert & Lacke, 2010), which is based 
on the famous Pareto-observation (Marshall, 2007). The most important lack of this method is 
that it cannot take more than one criterion into account at the same time. 

In recent years a number of multi-criteria decision models have been developed in order 
to compensate this defect of the traditional approach using methodologies such as the weighted 
linear optimization (Ramanthan, 2006), the genetic algorithm (Guvemir & Erel, 1998), the 
analytic hierarchy process (AHP) (Partovi & Burton, 1993), the fuzzy set theory (Chu, Liang & 
Liao, 2008), the fuzzy-AHP (Cakir & Canbolat, 2008), etc. 

On one hand these models can serve as management tools for solving complex decision 


situations (Vetschera, Chen, Hipel, & Kilgour, 2010), on the other hand the expansion of the 
number of dimensions might cause confusion in interpretation of the results. Not understanding 
or misunderstanding the decision context is one of the main problems (Ma, 2011) that can lead 
to mistakes. In general, the visualization of data can relieve the comprehension of the given 
problem and its environment surrounding it. 

Previously Tufte (1983) emphasized the importance of graphical representation analyzing 
the spread of cholera epidemic. He stated that the graphical analysis was more efficient than 
the calculation. Meyer (1991) examined the role of visual data in the concept of organizational 
research and suggested research questions, where visual data can be more meaningful to verbal 
one. Condon, Golden & Wasil (2003) proposed a methodology based on Sammon map for 
AHP. They analyzed the group decision procedure and relative judgment of the members. They 
claimed that using graphical assistant can raise the objectivity of collective decision making. 
Adler & Raveh (2008) introduced a graphical method for data envelopment analysis (DEA) using 
Co-plot, which can help to identify efficient units and outliers. Ma (2011) made up a screening 
methodology supported by graphical visualization based on the concept of multidimensional 
scaling (MDS). 

On the field of inventory management concerning the analyses of SKUs one of the main 
problems 1s that the large number of units can lead to impracticable data set that can make the 
employ of the given model for decision makers (DMs) more complicated and can hinder the 
appropriate interpretation of the results. In this sense the Case-based distance model developed 
by Chen, Kilgour & Hipel (2008) excels from other models considering the structure of the 
model, because it involves an untapped opportunity in the view of interpretation of results. 
Namely the process of determining the parameters of surfaces separating each group of SKUs 
is built in the model and visualizing these surfaces can help to discover the magnitude of groups 
and the relationships between them. 

The main purpose of the research is to give a graphical support to Multicriterial ABC 
Analysis based on Case-based distance model. Hence first of all the process of the initial dataset 
generation will be discussed. It is followed by an introduction into the visualization of the 
model, which is needed to explain the main findings of the research. A 2000-point illustrative 
example serves as demonstration of the plotting method step by step. 


Methodology of Research 


The main aim of the visualization is to demonstrate, how the Case-based distance model 
works and to present its results. The necessary dataset was generated on random basis in a 
controlled model. In this sense the control model means that the intervals were created in order 
to represent the most analised criteria like unit value, weight and lead time. On the other hand the 
three criteria are reasonable because of the representation of the model in three dimensions. 


Data Analysis 


The main goal of the Case-based distance model is to sort the stock keeping units into 


groups based on characteristics of the units as the set Q of criteria (Chen, Kilgour & Hipel, 


2008). Let Tbe the set of the units (4 )being analyzed, where A’ eT 2 ol |g =A,B,C 


based ong, €Q| j =1,2,...,m. 


The minimum (= ) and the maximum (™) values taken up in positions A and 
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A~ form an interval (Chen, Kilgour & Hipel, 2008). The Euclidean distances taken from these 


points determine the position of the given unit on criteriong ,, 1.e. the distance from the 
extreme values of the interval. 


Let the set Z T, | g = A,B,C bea subset of 7, that involves the units which 


repg = 


r 
repg 


ee (4 )s a (- a )s ee (4° ) (Chen, Kilgour 


principally represent the set T, and Z,.,, €T, | g = A,B,C is one of these units, where 


& Hipel, 2008). 


The Case-based distance model is constructed for n dimensions, which adopted on three 
criteria becomes capable to be plotted as a three dimensional system of co-ordinates defined by 


the three elements of Q, where respectively radius vector q, represents thex, g, the y and 


q, the z axis. In the system of co-ordinates the radius vectors can be 


given by the co-ordinates of their end, so if the co-ordinates of a P point are q, ; q5 


and q; , then the radius vector heading towards the P point is q, (4; 39539; )- 


In this sense we find that the coordinates of radius vector g__ belonging to oe el 
=y 


on criterion q ; : 


_ If the upper limit of the interval is the point of comparison: 
g,,€""* (4 )-c, Ging ) 05 0) 
g,,0; c2" (4 res ie ) 0) 
g,, 0; 0; 4" es Ein) 


- If the lower limit of the interval is the point of comparison: 
9, €: Erne Fer" (4° ) 0; 0): 
7,,0; € Eline Fes" (4 ) 0}: 
7,05 0; ¢3 Ein Fe 4) 
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After launching the normalization factor d*""” the distance between the lower and upper 
bounds of the intervals changes to unit, 1.e. (Chen, Kilgour & Hipel, 2008) 


1. Equation 


ox (greg) ke be 
dj (4 ,A je AE SS a1] f= 1,23 
dj 
Since the distance between the lower and upper bound of interval is equal to the absolute 


value of radius vector on criterion j, therefore the Euclidean distances taken from the points of 
comparison can be identified as the length of the radius vectors belonging to each unit. So the 


distance of Ze ef , from the upper limit of the interval on the second criterion is: 


is (4 2, pM ao) WG Grebo D 
2 »“re _ max 7 7 
PS d, (("" ( io (4* bem (4) o)) 


0° +" be, Ein. ) +o) EPO re Gow 


a max = max 
2 > 


where 


2. Equation 


asm = em et) =m bem e) = Ger Oeben eer" Do) = |) =a" 


Based on Chen, Kilgour & Hipel (2008) the weighted aggregation formula of distances are: 





3. Equation 


Gn, Coe )2€ oe) 
D* Cre )= DAG shee = D0) 4s Brae J = 8 a 
JE JE J ) 


4. Equation 


(1 6 inw )»Cin =n I) 
D Cree =D” nari L5G Y= 0) we 
je Je J 


respectively. 
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Note, that depending on the point of comparison we get two different sorting problems. 
Based on this fact the visualization has to be separated likewise: 


- If cn (4 ) is the point of comparison, then we talk about minimum transformation 
in the minimum (distorted) space; 


- Ife (4° ) is the point of comparison, then we talk about maximum transformation 


in the maximum (distorted) space. 


In this sense the attribute “distortion” means that the values created by the value function 


d j (cr : ) from the original values are linear and nonlinear distorted and this value transformation 
affects the sets and the shape of their separating surfaces, too. 

The main aim 1s to plot the surfaces (ellipsoids) bounding the sets in the original space. To 
achieve this goal the linear and nonlinear transformations together with their effects have been 
discovered and are introduced it in the following section. The softwares, Lingo and MATLAB 
were used to calculate the necessary model parameters and to plot the surfaces. 


Linear and nonlinear transformations and their effects on the visualization 


Space 0: Let P'e€ R° be the point belonging to a SKU with the co-ordinates 
P' = i pi,pi) where 7 = 1,2,...,D and 


1 

m, Sp, <M,, 
i 

m, Sp, <M,, 


l 
m, <p, <M,. 


where (m, ,m,,m, ) and (M »M,,M, ) are the minimum and maximum values of the original 


Co = ordinates, ie. Cc. (cm "Tc, 3c. in) and 
max max , max , max 
C 6 Ce iG, ) 
Space 1: Linear transformation:, / = (e, Sas: ) = T(P), 

where 

_ Pym _ Pz —™, _ P3 —M, 

1M 2 3 
1 —M, 5 —mM, 3 —M, 


i.e. in the case of i = 1,2,...,D 
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5. Equation 187 





i i i 

i; Pm ij Py —M, ;  P;—M, 

é) ~~ ) é, > p) oe ~~ 5) 
M,-™m, M,-m, M,-m, 


consequently the linear distorted co-ordinates are: 
ae a en > oe i 
i = (c' ei,ei)= T(P ), 


As aresult we get a unit cube placed in ,,o” origin. 


Space 2: Minimum world (nonlinear transformation): & = (é, oc. ) = R(E ) 


Let 
6. Equation yp (pi-m,) 


(Mf, —m,) 

5 e 

ea (eif - (Pia ima) 

(4, -m,)° 
a (ep = Pion) 

(4f, -m;) 
so in the case of 7 = 1,2,..., D oe =(é',é,€!)=3(tT(P’)). 
Note, that for 7=1,2,...,D O0<& <1,0<& <1,0<& <1, 


that is the transformed points ¢’ — q™ (2’ ) are in the same unit cube. 
j “repg 


Space 3: Maximum world (nonlinear transformation): 77 = (77, pedals ) == H(E ) 


7. Equation 
n= (nei) = Mine 
nm; =(I-2:) = ae 
m =(1-s3) = a 


so in the case of i= 1,2,...,.D  7' = (n' nimi) = H(T(P' )). 


Note, that in the case of 7 = 1,2,..., D 0<7 ae 0<7; an 0<7; <1, 
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that is the transformed 7 =a" (- r ) j — y are also in the unit cube. 
j repg 


It is easy to state that the relation between Space 2 and 3 1s: 


8. Equation 


Jé, tom =1, Jé +7, =1, JE, +n, =1. 


Depending on which space we are analyzing the position of the points P ' the parameters 
searched are the followings (Chen, Kilgour & Hipel, 2008): 


Space 2: &) =(E),6),8:)==(1(P') 
- weight vector: w = (w, WW; ); 
- radi: R,andR.. 
These parameters are given by the optimization model A MCABC expressed by the 
introduced transformation (& ) based on (Chen, Kilgour & Hipel, 2008): 
9. Equation 


min ERR=S'(ar Fh + Sark + (BF bE Sey, 
r=) r= ; 
where the conditions are: 


9b. Equation 


3 
—~ r r a Vee re . 
Sw, “Gy He SHG TH U2 ia5Me 


j=l 


9c. Equation 
3 
= r r -_ er me : 
yw, 6, PA, Si, | PH Lean, 
j=l 


9d. Equation 
3 
yw, G4 py ie | 2 ele 
j=l 


9e. Equation 
3 
yw, i Pg lg | FS arg 
j=! 
where 


Of. Equation 
0 Oe al lisa 0 ay | OM sae» Gn CO 
=e, 20. 02, 21 


-l<a, <0, 0<£f, <1, 
w. 20. dw, =1. 


jeEQ 
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Space 3: ii i 189 
P 1] = (yi .ni,ni)= H(T( ) 
- weight vector: w” = (w; ,W3 W3 } 
- radii: Rj and R,. 
Parameters are given in the case of A’ MCABC by the following model 
expressed by the introduced transformation (H _) based on (Chen, Kilgour & Hipel, 2008): 


min HRR = Ye) Ses) + (B,) Fe S(4e) ; 





10.Equation 


where the conditions are: 


10b. equation 


S] 














+ r r +]. ran. 
Sw POG SK, |\FHl2 cof, 
j=l 
10c. Equation 3 
+ r r +1. ran. 
Sw ‘Wh, +Q, SK, |r =152)...5n, 
j=! 
10d. Equation ; 
yw, oat Cie SOME cee i re 
j=l 
10e. Equation 
3 
+ r r + | ae 
yw Po SAP = lois 
j=l 
where 
10f. Equation 


Ven, a1 02k, a1) A, <i, 
-l<a’.<0, 0<f% <1 
-lo@,=0, U2 7, <1, 

w >0; yw el, 


jeQ 
Plotting method of linear and nonlinear separating surfaces 


The nonlinear transformation implied by the value function d j cr.) on 


Zo = Pong |i =r transformed the ellipsoids separating the set g into planes. In this way the 


planes separating the sets in the distorted spaces can be given by the following formula: 


Space 2: gi = (<',€:,é!)= =(T(P')) 
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11. Equation l 
é = —(R, - wg, -w26:), 
whereg=B ,C . 
Space 3: n' =(ni.n).ni)= H(T(P’)) 
12.Equation 





1 
+ + + 
mn (R —W, 7), —W, n,). 
ay 
"3 


a 
where g = A*,B”. 


Note that due the transformations above the assumed relation between Space 2 and 3 
cannot be assessed in distorted spaces. The DMs have to return back to the original space in 
order to determine the interaction of the two classifications in the same space. 

To get the three groups of SKUs and the ellipsoids separating them in the original space, 
executing the inverse transformation is needed. Since &' = (z',€:,€/)= S(T P ‘)) 


and n' = (n ; 1; : 1; ) = H(T(P )) are monotone transformations of P’ we can determine 
the ellipsoids separating the sets in the following way: 


Spee os (e221) =2{r(P')-> =(0(P')" =(0h ps pi)=P 


The equation of planes dividing the sets in Space 2 1s: 


a =e —W S, -w3é, ), 


W, 
which can be reformed using the 6. equation into following formula: 


GamY 1g, ng Com ym) 


(M,-m,) w; ° (M,—m,y (1, -m,) 


expressing P, can we find the equation of the ellipsoids: 






13.Equation 
P3 =M, 7+ Bary ay, amy) 
7 (M, me y 
where g=B ,C . 
Space 3: 


n’ =(n1.73.n5)=H(r(P')) > H(T(P'))" = (v1, ps. ps) =P 
The equation of planes dividing the sets in Space 3 is: 


| + + + 
i]5.= (ie WwW, 7), -win), 
Ww; 
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which can be reformed using the 7. equation into following formula: 


Hab) = 1g: -M,-py _,. M4, tp} 


(MV, —m, J W; “(y, —m,) (uv, —m,) 


expressing P, can we find the equation of the ellipsoids: 





14.Equation 





pe Ws =P) —p,J = MM, Pr) 


"=m wey 


where g = A”, B”. 


Note, that in this case the centre of ellipsoids is the upper limit of the available values, 
1.e. the maximum value. 


Since the position of the points P' representing the units in nm dimensional space 1s 
determined by the same criteria, therefore the sets are able to be plotted in the same space, 
where their sections can be identified as the nine sets before the reclassification. 


The linear and nonlinear transformations employed by the value function d , iG: — ) kept 
the relations between sets, that is why applying the 14. equation the reformed shape 
the planes separating the sets in Space 2 is 


15.Equation 


é,=1-|1- 





4/4 surface (,,eggs”), whereg = B,C. 
Respectively the reformed shape of the planes from the Space 3 can be given in Space 2 
by the following formula: 


16.Equation 





This surface is 4/4 surface (,,eggs”) likewise, where g = A”, B™. 
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Results of Research 





The results show that limitation of the 7 dimensional model into three dimensions made 
it possible to visualize the model in operation. The mathematical background of plotting method 
with the requisite formulas has been briefly introduced above. The following illustrative example 
supports to demonstrate the main findings of research concerning the visualization of Case- 
based distance model developed by Chen, Kilgour & Hipel (2008) step-by-step in figures. 

The illustrative example consists of 2000 units. Each unit has three parameters represented 
on the three axes of a three dimensional coordinate-system, unit value (x), ranges from 1,002 to 
17,385; weight (y) ranges from 9, 59 to 99,111 and lead time (z) ranging from 0,014 to 4,999 
assuming that the lead time is a continuos variable, which indicates the average of the regular 
lead times and in this sense it is not needed to be discrete. The intervals of unit value and weight 
were generated based on normal distribution, in the case of lead time a shorter interval was 
determined based on uniform distribution. 

Figure 1 illustrates the 2000-point dataset in three dimensional coordinate-system. The 
distribution of the unit value consists of three normal distributions with different parameters. 
Hence it is separated into three parts in order to make the movement of the dataset in the model 
to be traceable. 





Figure 1: Visualization of points P’ = (p! »P5>Ps ) i = 1,2,...,2000 in Space 0 
(original co-ordinates). 


Applying the parameters of the results of A -MCABC and A” MCABC each group 
can be represented in three dimensional 
coordinate-system. The model parameters are the followings: 


Space 2 (Minimum world): E' = (é' ,€, é!)= S(T(P' ) 
- weight vector: w” = (0,1840; 0,7241; 0,0919); 
- radii: R; =0, 2956 and Ro =0, 1508. 


Space 3 (Maximum world): 77! = (77 : 1; : 7; ) = H(T(P' ) 
- weight vector: w* = (0,1990; 0,7420; 0,0590): 
- radii: Rj =0, 2348 and R= 0, 4119. 
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Note, that both in the minimum and maximum worlds the centre of the comparison is_ | 193 
the “o” origin, but the order of sets is converse (see Figure 2). The converse order of these sets 
has been caused by the effect of linear and nonlinear transformations (introduced briefly 1n 6-7 
equations). 








Figure 2: Visualization of groups in Space 2 and 3 (marking: Black — group A, 
Grey — group B, Light grey — group C). 


Substituting the coordinates and the parameters above into equations 11 and 12 makes 
it possible to plot the surfaces separating each group from the others (see Figure 3). In both 
distorted spaces the separating surfaces quested for are planes. 





Figure 3: Visualization of groups and planes separating them in Space 2 and 
Space 3. 


The A, B and C groups can be plotted in the original spaces after executing the inverse 
1 ; ' 
transformation indicated as a(r(P’ y =a ry ope oo )= P' and 


H(t’) = (pi, De. pik P'. The bounding surfaces belonging to each group can be 
represented by substituting these reformed coordinates with changeless model parameters into 
equations 13-14 (see Figure 4). These surfaces in the original space are ellipsoids — nonlinear 
surfaces in accordance with the initial assumptions. 
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Figure 4: Visualization of groups of SKUs and ellipsoids bounding them based on 
the minimum and maximum world using the co-ordinates of 


x(t (P'))' = (pi, pi, pi) = Pi andH (T(P'))" = (pi, pi, pi)= Pi 


Due to the inverse transformation the results of the two classifications can be represented 
in the same original space (see Figure 5). 








Figure 5: Visualization of classification based on minimum and maximum worlds 
in one space (Space 0) — bottom-view and side-view. 


Since the implied linear and nonlinear transformations are monotonic, that is why the 
representation of the results of classifications keeps the relations between the groups in distorted 
spaces as well. The minimum and maximum worlds can be completed by using equations 15-16 
respectively (see Figure 6). 
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Figure 6: The reformed shape of the planes from Space 2 in Space 3 and from 
Space 3 in Space 2. 


Figure 6 shows that the shape of the surfaces in one world is oval (“egg’’) in the other 
world intercepting the planes. 


Discussion 


The adage “A picture 1s worth a thousand words” refers to the idea that a complex problem 
can be modeled in the way that managers can overview easily. The models like presented above 
can support every stages of decision making. Based on this managers will be able to make faster 
the process to get to the right decision. It increases the responsiveness of organizations. 

Visualization is extremly important in the case of Multicriterial decisions, when more 
factors are supposed to be taken into consideration at the same time. The plotting method 
proposed in this paper is a managerial tool to demonstrate, how the magnitude of groups varies 
caused by changing in the selected criteria. Knowing the linear and nonlinear transformations 
and the formulas of each surfaces the effect of criteria selection to the classification or the effect 
of switching from one criterion to another keeping the other criteria stable can be studied, which 
is one of the main contribution of the plotting method. 

In this sense the visualization can support not only the DMs but also the developers of 
the model to discover the opportunities in the model and to avoid the occurring problems. The 
visualization highlighted, what happens to the original data during the model that can make the 
explanation of results much easier. 

Discovering the movement and the distortion of the original dataset in each spaces could 
lead to the observation of the real relation between the two classification approaches (minimum 
and maximum worlds). Using formulas (equation 13-14) gained from the model allowed to 
produce Figure (4) which illustrates the position of the categories in the way as original model 
developers (Chen, Kilgour & Hipel, 2008) assumed. Furthermore, only the representation of 
the surfaces (ellipsoids) in the original space could help to determine the relationship firstly 
between the distorted spaces and secondly between the groups. 

On the other hand the visualization of the groups and surfaces in a three dimensional 
coordinate-system revealed that the rule of reclassification proposed by Flores & Whybark 
(1987) and used as a basis of reclassification in Case-based distance model (Chen, Kilgour 
& Hipel, 2008) is not unambiguos. Plotting the categories of SKUs from the two different 
approaches in one space still does not result in the final classification. The method of regrouping 
of the nine groups has to be reconsidered and can form a subject of further research. 
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Conclusions 


In recent years a number of models for Multi-criteria ABC analysis have been developed 
in order to make it possible to take more than one criterion into account at the same time. 
But in general, in the case of representing a model the introduction of the algorithm is more 
preferable than the visualization of data. However the visualization of results can relieve the 
comprehension of the given problem and the way how the model works even in relation to 
sorting problems based on Case-based distance model. 

The difficulty with interpretation of the results of classification is that the original 
data have been transformed by the value function employed in the model which affected the 
position of the groups and the shape of the surfaces bounding them. During the analysis of the 
mathematical background we contribute to the Case-based distance model with a graphical 
extension which is a proposal to the representation of the final classification of the groups of 
SKUs in one space. 
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